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1. Introduction and Motivation 

Redundant wavelet systems have been proved to be quit useful in many applica- 
tions, for examples, signal denoising, image processing, numerical algorithm, and so 
on. As a redundant system, it can possess many desirable properties such as symmetry, 
short support, high vanishing moments, and so on, simultaneously. In this paper, we 
are interested in the construction of tight wavelet frames with such desirable properties. 
Other than redundancy, computational efficiency is also desired in applications, which 
means a wavelet system should have as few generators as possible. It is well known 
that the construction of a wavelet system from a refinable function can be formulated 
as a matrix extension problem. In order to have fewer generators for a tight wavelet 
frame, we shall employ the algorithm for matrix extension with symmetry. 

Let us first introduce some notation and definitions. We say that d is a dilation 
factor if d > 1 is an integer. Throughout this paper, d is fixed and denotes a dilation 
factor. We say that : R i-> C is a d- refinable function if 

cf>^dJ]ao(md--k), (1.1) 

teZ 

where ao : Z i-> C is a finitely supported sequence on Z, called the low-pass filter 
for (p. The symbol of aq is given by ao(z) '■— Y^tez cioik)z!'. The coefficient support of 
ao is defined by coeffsupp(ao) := lin,n], where aoim)aoin) + and oqW - for all 
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k i [m, n]. In frequency domain, the refinement equation in ( II. Il l can be rewritten as 

^(d^)^ao(ma ^eR, (1.2) 
where oq is the Fourier series of gq given by 

The Fourier transform / of / e -Li(M) is defined to be /(^) = J^f(t)e^"^dt and can be 
extended to square integrable functions and tempered distributions. A low-pass filter 
ao is orthogonal if Z%o \ao{^ + 2^ ;/d)l^ = 1 for all ^ e K. 

Usually, a wavelet system is generated by some wavelet function ij/^ ,{ - 1, . . . , L 
from a d-refinable function as follows: 

-//(d^) = flK^)0(a ^ = (1.4) 

where each : Z i-> C is a finitely supported sequence on Z, called the high-pass 
filter fox ,€ = 1,...,L. 

We says that {i/r' , . . . , i/^^) generates a d-wave let frame in -L2(K) if {iA^^ := d^^^ifr^id^- 
—k) : j,k € 'Z,^ - I, . . . , L} is a frame in L2(R), that is, there exists two positive 
constants Ci , C2 such that 

L 

ciii/ii?,(R) < 2 Z Z i<^' '^i^>i' ^ c^ii/iiLk) ^/ ^ '^^(R), (1.5) 

f=l jeZ ieZ 

where \{f, tA^^)P = (.A '/';it)^^yi' /) ^'^'^ ') '^^e inner product in L2(M) defined by 

Jr 

If Ci = C2 = 1 in dl.Sb . we say that {1^', . . . ,1/'^) generates a f/g/jf d-wavelet frame in 
L2(K). 

An important property of a wavelet system is its order of vanishing moments. We 
say that {i^', . . . , i/'^) has vanishing moments of order n if 

J t''if/^(t)dt = A; = 0,...,n- l,Vf = 1,...,L. (1.6) 

(II.6I 1 is equivalent to saying that ^(f/'^iO) = for all ^ = 0, . . . , n - 1 and i = 1 , . . . , L. 

Let be a compactly supported d-refinable function in L2(]R) associated with a 
low -pass filter ao such that ^(0) - 1. Suppose there exist high-pass filters ai, . . . , 
such that 

L 

Y^aeaei- +2nk/d) = 6k, = 0, . . . , d - 1, (1.7) 

f=0 

where 6 is the Dirac sequence such that 6{Q) = 1 and d{k) - for all k Q. Define 
, { - 1 , . . . , L as in ( 11.41 ). Then {1/'' , . . . , i/r^) generates a tight d-wavelet frame, see 
Ullfioll . Moreover, if the low-pass filter oq satisfies 

1 - |5J,(^)|2 = (9(|#«), (1.8) 
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which means 1 - Iflo©!^ has zero of order 2n near the origin, then {(/r' , . . . , iff^} given 
by ( 11.41 ) has vanishing moments of order n, see |4]. 

Now, we are ready to introduce our main resuh of this paper on the construction 
of symmetric complex tight d-wavelet frame. Throughout this paper, we shall denote 
Pm.niy) a polynomial of degree « - 1 as follows: 

«-i 

Pm,n(y) = ^ 
j=0 



d-1 



z n 

./i+--+;rt-i = / k=l 



- 1 + jk 
jk 



sin(^;r/d) 



-Ik 



(1.9) 



Theorem 1. Let d be a dilation factor Let m,n e ^ be positive integers such that 
2n — 1 ^ m. Let Pm,n(y) be the polynomial defined in (11.9b . Then 

Pman-i(y)>0 VyeM. (1.10) 

Let zi,Zi, ■ ■ ■ ,Zn-\,Zn-\ €C\Wbe all the complex roots ofP,„2n-i- Then 

Pm,2n-l(y) ^ \Q,n,n(y)\\ (LH) 

where Qm,niy) - c(y-Zi) ■ ■ ■ (y-Zn-i) withe - (-l)""'(zi ■ ■ -Zn-iT^- Define a low-pass 
filter flo by 

flo(^):=e'L— J^l _ 1 e,„„(sin2(^/2)), (1.12) 



where [■} is the floor operation. Then 

m(d - I) 

= e'(i-«)f55(^) with e = OT(d - 1) - 2L— — -J, (1.13) 

and 

mid — 1) m(^ci — 1) 
coefrsupp(flo) = [-L ^ J - n + 1, L J + n-l+s]. 

Let (p be the standard d-refinable function associated with the low-pass fllter aQ, that is, 
0(^) Y\J=i '^)(d ■'^)- Then (p is a compactly supported d-refinable function in L2(K) 
with symmetry <^(grf ~ ) — 4'- Moreover, one can construct high-pass filters 
L € {d — l,d, d + 1), with symmetry by Algorithm^ such that ( 11.71 ) holds. Define 
as in ( |1.4| l. Then [ip^ , . . . generates a tight d-wavelet frame in L2(M) and 
{(/'', . . . has vanishing moments of order 2n — 1. Furthermore, |coeffsupp(flf)| ^ 
|coeffsupp(flo)| and (/'^(jrf + '7^ - = B^i]/^ for some rf e {1,0, -1) and e {-1, 1}, 
( — 1 , . . . , L. In particular, for m = 2n — 1, L = d-1 and {ip^ , . . . , (A'^"' ) generates an 
orthonormal wavelet basis with symmetry for L2(K). 

Theorem[T]is summarized from the following sections. In section 2, we shall in- 
spect the important properties of P,„,„(3') and show that for any dilation factor d and 
m, n e N, we can construct a symmetric complex d-refinable function from Pm,2n-i{y), 
called complex d-refinable pseudo spline of type I with order (m,2n - 1). Moreover, 
we prove that the shifts of all such complex pseudo splines are linearly independent. 
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In section 3, we introduce the general problem of matrix extension with symmetry and 
present a step-by-step algorithm for the matrix extension with symmetry. In section 4, 
we shall apply our matrix extension algorithm and show that we can construct symmet- 
ric complex tight d-wavelet frames from symmetric complex pseudo spline with only 
d wavelet generators. Moreover, we shall provide several examples to demonstrate our 
results and algorithms. Conclusions will be given in last section. 



2. Symmetric Complex Reflnable Functions from Pseudo Splines 

In this section, we shall restrict our attention to the construction of d-refinable 
functions from which wavelet system with high vanishing moments can be derived. In 
other words, we need to design low-pass filters a such that a satisfies ( 11.81 ) for some 
n e N. 

To guarantee that the d-refinable function associated with a has certain regularity, 
usually the low-pass filter a satisfies the sum rules of order m for some m e N: 



(1 + e-'^ + ---+e 



-/(d-l)^\m 



That is, a is of the form: 



1 + e-'f + ■ ■ ■ + e-'(d-i)f 



(2.1) 



(2.2) 



for some 2;r-periodic trigonometric polynomial L with L(0) - 1 . For L(^) = 1 . a(^) is 
the low-pass filter for B-spUne of order m: B^(£) = e'L'«/2Jf(i _ e-'¥)'«/(/^)'«. 
Define a function h by 



d-i , 
k=i ^ 



ye 



We have 



sin {kn/d)l 
1 + --- + e-'<^-i)f|2 sin2(d^/2) 



d2 sin2(^/2) 



and 



where 



Ky)- 



6-1 



nz 



^0 sin^Hkn/d) j 



d-I 



z n 

jl+-+jd-l=j k=\ 



m-l+ jk 
Jk 



J] CmjyK \y\ < sin^TT/d), 

i=0 



sin(/t7r/d)"^-", e N 



(2.3) 



(2.4) 



(2.5) 



(2.6) 



Note that Pm,n(y) - Tj'j=o '^mjy-'- Consequently, we have the following result. 

Lemma 2. Lettn,n e N be such that n ^ m; let f ,„ „ and h be polynomials defined as in 
(11.9b and (12.31 1, respectively. Then f ,„ „(sin^(^/2)) is the unique positive trigonometric 
polynomial of minimal degree such that 



1 -/i(sin2(^/2))"P„,„(sin2(^/2)) = 0(1^1^"). 



(2.7) 
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Proof. Suppose there exist trigonometric polynomial g(^) such that 



1 - h(smH^/2)rg(0 



Observing that the first (n - 1 )f/z-degree Taylor polynomial of h(y) 
we have 



' aty = Ois Pm,n(y)' 



g(0 = h(sm\^/2)y 



"[1 + 0(\^\^")] 



P,„,„(sin2(^/2)) + (9(|sin2(^/2)|"). 



Moreover, it is easily seen that the coefficients of Pm^„{y) are all positive. Consequently, 
f ,„_„(sin^(^/2)) is the unique positive trigonometric polynomial of minimal degree such 
that dZTb holds. ■ 

For m,n e N such that n ^ m, let TJai^) := /!(sin^(^/2)'"P,„,„(sin^(^/2)). Then the 
d-refinable function //0 associated with //« by dl.ll ) is called the pseudo spline of type 
II. By Lemma|2l using Riesz Lemma, one can derive a low-pass filter /o from //o such 
that = //fl(^). The d-refinable function /0 associated with such /o by (11.11 ) is 

referred as real pseudo spline of type I. Interesting readers can refer to [? ? ? ? ] for 
more details on this subject for the special case d = 2. 

Note that /« satisfies dl.SI ). One can construct high-pass filters b^,...,b^ from 
/fl such that ( 11.71 ) holds. Then t//^ , . . . , t//^ defined by ( 11.41 ) are real-valued functions, 
{i/'', . . . , i/r^) has vanishing moment of order n and generates a tight d-wavelet frame. 
However, {i/r', . . . , t//^} is not necessary symmetric since the low -pass filter ja from //« 
via Riesz lemma might not posses any symmetry pattern. In what follows, we shall 
consider complex-valued wavelet generators and show that we can achieve symmetry 
for any odd integer n e N. 

The next lemma is needed later, which generalizes [? , Lemma 2.1], 

Lemma 3. Let P,„,„(y) be defined as in il.9\ . Then 

d-1 



{=1 



(1 



P,n,n(y) 



(2.8) 



where 



>l+'-+>d-l='l-l 

Proof. By (11.91 ), we have 



m + ]( 
jc 



sin(-)-20Hi 
d 



d-1 

■ n 

k=lM[ 



P'm,n(y) 



n-2 

Z 



d-1 



n-2 



d-1 



J=0 



E n 

<-7d-i=;+i *=i 

Z Z 

■■+;d- 

Z 



- 1 + jk 
Jk 

d-1 



sin(5)-2>' 
d 



- 1 + A 

Jk 



(; + i)y 



sini — y^-"' 
d 



-n 

*;=1 



- 1 + Jk 
Jk 



sin(-)-2>' 
d 



d-1 n-2 

f=l y=o 
d-1 

=: m 

>=i 



;i+--+yd-i=./ 



m + j( 

Je 



sin(-) 



20V+1) 



d-1 

n 

k=\MC 



m-\+ jk 

Jk 



sin(— )"^^' 
d 
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where 



n-2 



QKy) X X 



m + je\ . 
I sin' 

jc 



d-I 



(T>""'*" n 



in - 1 + jk\^._,kji^_2j^ 



Jk 



sin(-)- 



Note that 

y 



sin^(f) 

n-2 

=z z 

>=o L>i+--+id-i=j 
«-i r 



je 



d-l 

sin(--)-20V+2) ]-[ 

l(:=I,l(:^f 



Jk 



sin(-)- 



z z r; 'r>'"<f n 

,■-1 ;.4-...j-ij . = ; V .11 I " i— I i.. 



d-l 



k=\MC 



j=i L)i+--+;d-i=./ 

Using ("'+/^) - ('"T/_y^) = ('""];^'^), it is easy to deduce that 



m -I + jk 
jk 



sin(— 
d 



y. 



I sin2(^)j 



sin^(^) 



where 



Z 



;i+'-+;d-i=«-i 
Consequently, ( I2.8l l holds 



sin(— r^^''^' 
jc I d 



d-l 

' n 

k=lMl 



^Msin( — 
Jk I d 



By the above lemma, we have the following result regarding the positiveness of 
Pin,n{y)^ which generalizes [? , Theorem 5] and [? , Theorem 2.4] and proves ( 11.101 ). 

Theorem 4. Let m,n e N be such that n ^ m. Then P,n,n(y) > 0/or a// y e R if and 
only ifn is an odd number. 

Proof. It is easily seen that Pm,n(y) > for all y due to the positiveness of its 
coefficients and limy^_co Pin,n(y) - when n is even. Hence n must be an odd 
number if Pm,n(y) > for all y e M. 

On the other hand, suppose n is an odd number Let yo < be any stationary point 
of Pm.niy)- Then by Lemma[3] we have 



d-I 



Pm,n(yo) ( „_i 
Sin (f ) 



which implies that Pm,n(yo) > 0. Consequently, Pm,n(y) > for all y < 0. This 
completes our proof. ■ 
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Now, by Pm,2n-i(y) > for all y e R and 2n - 1 ^ m, P,n,2n-i(y) can only have 
complex roots. Hence, we must have 

n-l 

Pm.2n-i(y) = Co Y\(y ' Zj){y - Tj), zi,T[, ■ ■ .,z„-uzii^ e C \ R. 
;=1 

Define Q,„j,(y) = c Hpi ~ Zj) as in dl.lll i and flo(^) as in ( 11.121 ). Then it is easy to 
check that the symmetry pattern of oq satisfies ( |1.13t . We shall refer the d-refinable 
function 0„,_„ associated with the low -pass filter ao as complex pseudo spline of type I. 

Now, we have the following result which shall play an important role in our con- 
struction of tight d- wavelet frame in section 3. 

Corollary 5. Let m, n e N be such that 2n - 1 ^ m and ao be defined as in (11.12b . 
Then 

d-1 

1 - \^'>(^ + 2'^;/cl)|' = \boim\ (2.9) 

/or some In-periodic trigonometric function bo(^) with real coefficients. In particu- 
lar, bl"-''"(^) = and = C2„,2«-i[sin2(^/2)/d2]2«-i with C2„,2«-i being the 
coefficient given in i2.6i . 

Proof. We first show the 1 - Z%o \ao{^ + 27r;/d)p > for all ^ e M. Let yj 
sin2(^/2 + nj/d), j = 0, . . . , d - 1. Since = h(sm\^/2)y"P„,2n-dsm^(^/2)), we 

have 

d-1 d-1 
1 - 2 + 2nj/dt = 1 - h(yjrP,n,2n-l(yj) 

d-1 d~l m-1 

= 1 -J]hiyjTP^Ayj) + Yj^iyiT ^ 

;=0 j=0 k=2n~\ 

d-1 m-I 

;=0 k=2n-l 

The last equality follows from the fact that the low-pass filter a satisfying = 
h(yo)'"Pm,m(yo) is an orthogonal low-pass filter , which satisfies Yi'jZo h(yj)'"Pm,m(yj) = 
1 ([? ]). Now, by that 1 - X%o l«o(^ + 27r;/d)p ^ is of period 27r/d, ^ follows 
from Riesz Lernma^ 

Obviously, = since fly""' " is an orthogonal low-pass filter. And from 

above, noting that h(yj)yj - sin^(d^/2)/d^. we have 

d-1 d-1 

\bl"''-"-Um^ = C2„,2„-i ^/<3';)'"3'f-' = C2„,2„-i Y^h(yj)yjf"-'Myj) 
j=o j=o 

d-1 

= C2„,2n-dsinHd^/2)/d^f"-'Y>'(yj)Pli(yj) 

;=() 

= C2„,2„-i[sin2(d^/2)/d2]2"-i. 
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Next, we shall discuss the linear independence of the pseudo spline of type I (real 
and complex) and type II. Linear independence is an important issue in both approx- 
imation theory and wavelet analysis. The linear independence of the integer shifts of 
a refinable function is a necessary and sufficient condition for the existence of a com- 
pactly supported dual refinable function, see [? ]. For a compactly supported function 
(f) e -L2(]R), we say that the shifts of are linearly independent if 

2cOX--;) = implies c{j)^Q^je1. (2.10) 

We say that the shifts of (p are stable if there exists two positive constants C\ and C2 
such that 

Ci2|c(;)|2 ^ ||^c(;X--;)||2^(jj) ^ C2 ^ |c(;■)|^ (2.11) 

jeZ jel. jel. 

for all finitely supported sequences c : Z i-> C. 

It is known in [? ] that the shifts of a compacly supported function (p e L2{M.) is 
linearly independent (or stable) if and only if 

span{^(^ + Ink) : ;t e Z) ?t for all ^ e C (or for all ^ e R). (2.12) 

Here, the Fourier-Lapalce transform of <p is defined to be 

0(^)= 1 if,(t)e-'^'dt, ^eC. (2.13) 
Jr 

Consequently, if the shifts of a compactly supported function (p e LaiTS.) are linearly 
independent, then the shifts of (p must be stable in L2(IR)- 

The following lemma is needed later to prove the linear independence of the shifts 
of pseudo spline of type I or type II. 

Lemma 6. Let m, n e N be such that I < n < m. Let Cmj be the coefficient of P„,^n(y) 
defined in ( 12.61 ). Then 

2c„j-i < c„,j, = 1, . . . n - 1. (2.14) 
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Proof. First, it is easy to show that 2('" < ('J^/) for ; = 0, . . . , n - 2. Then 

2..,,,- =2 2 n("'";;>'"<T'"" 

d-l / , • x „ d-l 



< J_ V y + f ) sin(^)-OV.i) pj /'« - 1 + ;A 



2A 



2>j 



Due to = /0(^) ■ i4'{-^) for all ^ e C, the shifts of the pseudo spline of type 

II are linearly independent will implies the shifts of the pseudo spline of type I are 
linearly independent as well [? , Proposition 1.1]. Moreover, we have the following 
lemma, whose proof is also similar to [? , Lemma 2.1] 

Lemma 7. Let e LiiK) be a compactly supported refinable function associated with 
a low-pass filter a by (ILll l. Then, the shifts of 4> are linearly independent if and only if 

(1) the shifts of (p are stable; 

(2) the symbol a of a satisfies: {aize'^"^"^) : ; = 0, . . . , d - 1) 0/or all z e C. 

For a polynomial P{z) - cq + ciz + ■ ■ ■ + c„z" with real coefficients satisfying 
c„ > c„_i > • ■ • > Co > 0, it was shown in [? ] that all zeros of P(z) are containing 
inside the open unit disk {z € C : |z| < 1). 

Using the above facts, we can prove the following result. 

Theorem 8. The shifts of the pseudo spline of type I (real or complex) or type II are 
linearly independent. 

Proof. We only need to show that the shifts of the pseudo spline of type II are linearly 
independent. 

Fix m, n e N with 1 ^ « ^ m. Let the low-pass filter a be given as follows: 
a(0 = /z(sin2(^/2))'"P„,,„(sin2(^/2)). 
Let <p be the pseudo spline of type II associated with mask a by dl.ll l. 



9 



For n - m, (p associated with a is an interpolatory d-refinable function [? ], i.e., 
(pij) = 6{j), 7 e Z and 0(d^) = a(^)4>{^). In this case, the shifts of must be Hnearly 
independent. 

For n - 1, by that the shifts of the B-spHne of order 2m are stable and the symbol 
of a is a(z) = |1 + ■ ■ ■ + z''"'p"Vd^"', the fact that the shifts of (p associated with a are 
linearly independent follows by Lemma|7] 

Hence, without loss of generality, we can assume that \ < n < m. Since a(^) ^ 
/!(sin^(^/2))"' and the shifts of B-spine of order 2m are stable, by ( |2.12t . the shifts of 
(p are stable. We need to prove that the symbol of a satisfies item (2) of Lemma|2] By 



sin^(^/2) = with z = e"'^, we have 



a(z) = 





1 + z + 










d 






1 + z + 




+ z^-' 






d 






1 + z + 








d 



2m n-l 



2m n-l 

i=0 



4z 



2z 



Now, item (2) of Lemma [7] is equivalent to saying that there is no z e C such that 
P(^gi2;rj/d) = for all = 0, . . . , d - 1, where P(z) 
not, let zo = pe'" e C be such that P(zj) - for all ; 



2-1 j=o ^ '-I 



■ Suppose 
, d - 1, where Zj - pe'^' 



with 61; = 6' + 27r77dfor j 



0, . . .,d - 1. Note that 



2p 



2p 

2zi 



cos(0;) ^ 1 - cosiOj), = 0, . . . , d - 1 . 



^ 1 . This contradicts to the fact that all roots of 



Consequently, maxo^j^cj_i |- 
the polynomial Y!]Z{)'^ ''cm,iZ' are containing inside the unit disk {z e C : |z| < 1) 
due to its coefficients 2^Cm,j,j - 0, . . . ,n - 1 satisfying 2^^^'^Cm,j-\ < 2^^Cmj for all 
7 = l,...,n-l. We are done. ■ 



3. Matrix Extension with Symmetry 

In this section, we shall introduce the matrix extension problem, which plays an 
important role in our construction of tight d-wavelet frame. In order to state the matrix 
extension problem and our main results, let us introduce some notation and definitions 
first. 

Let p(z) = Yjkez Pkz'^^ z 6 C\{0) be a Laurent polynomial with complex coefficients 
pic & C for all k e Z. We say that p has symmetry if its coefficient sequence {pk}k£Z has 
symmetry; more precisely, there exist e € {-1, 1) and c € Z such that 

Pc-k^epk, V^eZ. (3.1) 
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If e = 1, then p is symmetric about the point c/2; if e = -1, then p is antisymmetric 
about the point c/2. Symmetry of a Laurent polynomial can be conveniently expressed 
using a symmetry operator S defined by 

Spiz) := z e C\{0). (3.2) 

When p is not identically zero, it is evident that ( 13.1b holds if and only if Sp{z) = 
ez'^ ■ For the zero polynomial, it is very natural that SQ can be assigned any symmetry 
pattern; that is, for every occurrence of SQ appearing in an identity in this paper, SQ is 
understood to take an appropriate choice of sz^ for some e e {-1,1) and c e Z so that 
the identity holds. If P is an r x 5 matrix of Laurent polynomials with symmetry, then 
we can apply the operator S to each entry of P, that is, SP is an r x s matrix such that 
[SP} j^k '■- S{{P]j,k), where [P];.* denotes the {], A;)-entry of the matrix P throughout the 
chapter. 

For two matrices P and Q of Laurent polynomials with symmetry, even though all 
the entries in P and Q have symmetry, their sum P + Q, difference P - Q, or product 
PQ, if well defined, generally may not have symmetry any more. This is one of the 
difficulties for matrix extension with symmetry. In order for P ± Q or PQ to possess 
some symmetry, the symmetry patterns of P and Q should be compatible. For example, 
if SP = SQ, that is, both P and Q have the same symmetry pattern, then indeed 
P ± Q has symmetry and S{P + Q) - SP - SQ. In the following, we discuss the 
compatibility of symmetry patterns of matrices of Laurent polynomials. For an r x s 
matrix P(z) - Yjkez Pkz'', we denote 

P*(z):= with Pl:^Tk, keZ, (3.3) 

keZ 

where f ^ denotes the transpose of the complex conjugate of the constant matrix Pk in 
C. We say that the symmetry ofP is compatible or P has compatible symmetry, if 

sp(z) = (seiy(z)S02(z), (3.4) 

for some 1 x r and I x s row vectors 0i and 02 of Laurent polynomials with symmetry. 
For an r X s matrix P and an 5 x f matrix Q of Laurent polynomials, we say that (P, Q) 
has mutually compatible symmetry if 

SPiz) ^ iS0inz)S0iz) and SQiz) ^ (S0nz)S02{z) (3.5) 

for some Ixr, Ixi, Ixf row vectors 0i , 0, 02 of Laurent polynomials with symmetry. 
If (P,Q) has mutually compatible symmetry as in ( 13.51 ), then it is easy to verify that 
their product PQ has compatible symmetry and in fact >S(PQ) = (S0iyS62- 

For a matrix of Laurent polynomials, another important property is the support of 
its coefficient sequence. For P = Ejtez Pkz!' such that Pk - for all k e Z\[m, n] with 
P,„ + and P„ + 0, we define its coefficient support to be coeffsupp(P) := [m, n] 
and the length of its coefficient support to be |coeffsupp(P)| :- n - m. In particular, 
we define coeffsupp(O) := 0, the empty set, and |coeffsupp(0)| :- -oo. Also, we use 
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coefF(P, k) := to denote the coefficient matrix (vector) Pk of z*^ in P. In this chapter, 
always denotes a general zero matrix whose size can be determined in the context. 

The Laurent polynomials that we shall consider in this paper have their coefficients 
in a subfield F of the complex field C. Let F denote a subfield of C such that F is closed 
under the operations of complex conjugate of F and square roots of positive numbers 
in F. In other words, the subfield F of C satisfies the following properties: 

xeF and ^Jye¥, Vx,yeF with y > Q. (3.6) 

Two particular examples of such subfields F are F = M (the field of real numbers) and 
F = C (the field of complex numbers). A nontrivial example other than R and C that 
satisfies i3.6i is the field of all algebraic number, i.e., the algebraic closure Q of the 
rational number Q. A subfield of R given by Q n R also satisfies ( 13.6b . 

Now, we introduce the general matrix extension problem with symmetry. Let r 
and s denote two positive integers such that 1 ^ r ^ Let P be an r x matrix of 
Laurent polynomials with coefficients in F such that P(z)P*(z) - Ir for all z e C\{0} 
and the symmetry of P is compatible, where Ir denotes the r x r identity matrix. The 
matrix extension problem with symmetry is to find an i x i square matrix P^. of Laurent 
polynomials with coefficients in F and with symmetry such that [1^, 0] Pg = P (that is, 
the submatrix of the first r rows of P^. is the given matrix P), the symmetry of P^ is 
compatible, and Pf.(z)PJ(z) = Is for all z e C\{0) (that is, P^ is paraunitary). Moreover, 
in many applications, it is often highly desirable that the coefficient support of P^ can 
be controlled by that of P in some way. 

In the context of wavelet analysis, matrix extension without symmetry has been 
discussed by Lawton, Lee and Shen in their interesting paper ^ and a simple algorithm 
has been proposed there to derive a desired matrix P^ from a given row vector P of 
Laurent polynomials without symmetry. In electronic engineering, an algorithm u sing 



the cascade structure for matrix extension without symmetry has been given in 1112 1 
for filter banks with perfect reconstruction property. The algorithms in [8, lij] mainly 
deal with the special case that P is a row vector without symmetry and the coefficient 
support of the derived matrix P^ indeed can be controlled by that of P. For F = R and 
r - \, matrix extension with symmetry has been considered in [9] . 

We study the general matrix extension problem with symmetry and have the fol- 
lowing results. 

Theorem 9. Let P be an r X s matrix of Laurent polynomials with coefficients in a 
subfield ¥ of C such that ( 13.6b holds. Then P(z)P*(z) = Ir for all z e C\{0) and 
the symmetry of P is compatible as in (13.4b . if and only if there exist s X s matrices 
P(), . . . , Py+i of Laurent polynomials with coefficients in F such that 

(1) P, := Py+i(z)Py(z) ■ ■ ■ Pi(z)Po(z) is paraunitary: P,(z)P:(z) = h; 

(2) P can be represented as a product of Po, . . . , Py+i." 

P(z) = [/„0]Py+i(z)Py(z)---Pi(z)P,)(z); (3.7) 

(3) Pj, 1 ^ /' ^ 7 are elementary: Pj(z)P*(z) = ant/ coefFsupp(Py) C [-1, 1]; 

(4) (Pj+i, Pj) has mutually compatible symmetry for all ^ y ^ J; 
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(5) Po - Ugg^ and Pj+\ - diag{\Jsei,h-r), where \^se^, ^soi ^'"^ products of a 
permutation matrix with a diagonal matrix of monomials, as defined in ( 13. 9I )." 

(6) J ^ max [|coeffsupp([P],„„)|/2], where \-~\ is the ceiling function. 

The proof of Theorem |9] is given in [? ], in which a general extension algorithm 
(r ^ 1) for the matrix extension with symmetry was also proposed. Since we are 
only concerning about the construction of tight d-wavelet frame in this paper, we shall 
provide an simple algorithm for the matrix extension with symmetry for the case r - \. 

Before we introduce our algorithm, we need some auxiliary results. 

For a 1 X n row vector f in F such that ||f|| + 0, where ||£|p := ff*, we define 
«f to be the number of nonzero entries in £ and := [0, . . . , 0, 1, 0, . . . , 0] to be the 
y'th unit coordinate row vector in M". Let be a permutation matrix such that fiif = 
[/i , . . . , /,„ , 0, . . . , 0] with /j ^£ for ./■ = 1 , . . . , nf . We define 

f ill. if«f = 1; 

Vf:={ , . (3.8) 

where Vf := £ - ^l|f||ei. Observing that ||vf||2 = 2||£||(||f|| - |/i|), we can verify that 
VfVf = /„ and fijfVf = ||f||ei. Let f/f := /SfVf. Then f/f is unitary and satisfies 
Ui = [§;\.F''] for some (n - 1) x n matrix F in F such that f f/f = [||f ||, 0, . . . , 0]. We 
also define f/f := /„ if f = and t/f := if f = 0. Here, f/f plays the role of reducing 
the number of nonzero entries in £. 

Let be a 1 x « row vector of Laurent polynomials with symmetry such that SO - 
[siz"^' , . . . , e„z''"] for some si, . . . , e„ € {-1, 1) and ci, . . . , c„ e Z. Then, the symmetry 
of any entry in the vector 6'diag(z"'^'^' . . . ,z^^'^"^^'^) belongs to {+1, Thus, there 

is a permutation matrix Eg to regroup these four types of symmetries together so that 

S(e\^se) = [l„,,-U,z-'l„„-z-il„J, (3.9) 

where \Jse '■- diag(z"^'^' . . . ,z''^'^"^^'^)Eg, 1„, denotes the 1 x m row vector [1, . . . , 1], 
and n4 are nonnegative integers uniquely determined by SO. Note that U^e do 

not increase the length of the coefficient support of 9. 

Let q be a 1 X 5 row vector of Laurent polynomial satisfying qq* = 1 and >Sq - 
[lij, -l,2,z"'ls,, -z^'ljj] for some nonnegative integers si, . . . , S4 such that s\ + S2 + 
S3 + S4 = s. Then, q must take the form in ( 13. 10b or ( 13.11b with f 1 as follows: 



q = [fi, -£2,91,-92]/' + [£3, -£4, 93, -94]/'^' + coeff(q,^)/ 

f=fl+2 

+ [£3, £4, 91, 92k''"' + [f 1, £2, 0, Ok'^ 



(2-2 



q = [0,0,£,,-£2]z'> + [9i,-92,£3,-£4k''^' + ^ coeff(q,^)z' 



c=e, +2 



(3.10) 



(3.11) 



+ [93, 94, £3, £4]/'"' + [91, 92, £1, £2]/'. 
If q takes the form in ( 13.1 lb . we further construct a permutation matrix Eq so that 
[gi, 92, £1, fil-Eq = [fi, £2, gi, 92] and define Uq,£ := Eqdiag(Is-s^,z'^Is,), where Sg is 
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the size of the row vector [gi, 92]- Then qUq takes the form in ( 13.101 ). For qZSg of form 
(13.101 ), we simply let Uq := /j. In this way, qo qUq always takes the form in ( 13.101 1 
with £1 0. 

Note that UqUq = Is and ||£i|| = ||£2|| if qoqg - 1. Now an s x paraunitary matrix 
Bq,, to reduce the coefficient support of qo as in ( 13.101 ) with fi from {{1,(2} to 
Ui + IJ2- 1] is given by: 





cFi 





gi(i + ^) 



92(1 -i) - 



1 


-fi(z-7) 



CF2 


-91(1-7) 



-92(1 + i) 



c 


^7^1(1+^) 

<.fi 







cGi 











-Jf2(l+Z) 








cg;g2 

cG2 



(3.12) 



where C£, := ||f Cg, := ||gi||, Cg^ := ||g2||, cq := coeff(qo, £1 + l)coeff(q;, -^2)/cf,, 

— T — if gi ^ 0; —J — if 92 i= 0; 

c otherwise, " otherwise, (3.13) 

c:=(44 +24 +2cl_ + \co\y\ 

and = t/f^, [g",G*] = f/g, are unitary constant extension matrices in F for 

vectors fj, gj in F, for j = 1,2, respectively. The operations for the emptyset are 
defined by ||0|| = 0, + A = A and • A = for any object A. 

Lemma 10. Let q be a I x s row vector of Laurent polynomial satisfying qq* — 1, 
|coeffsupp(q)| > 2, andSq - [1^,, -Ij^z^'lj,, -z"'l.sJ/or wme nonnegative integers 
si, . . . , S4 such that si + S2 + sj, + S4 — s. Let Bq :— UqBq(,U* with Uq and Bq„ being 
constructed as above. Then, >SBq = [1.5,,-1.5j,zlj3,-zljJ^[lj,,-lj,,z"'lj,,-z"'1.5j, 
coeff'supp(Bq) = [-1, 1], and coefFsupp(qBq) - {(\ + 1,^2 - 1]- That is, Bq has com- 
patible symmetry with coefficient support on [-1,1] and Bq reduces the length of the 
coefficient support o/q exactly by 2. Moreover, >S(qBq) — Sc\. 

Proof. Direct computation shows that coeff(qB, k) = Q for k e {t\ - \,l\,{2 + 1, ^2)- 
One can also show that coeff(qB, fc) ^tQforfce {^1 + 1,^2-1)- Hence coefFsupp(qBq) - 
1 + 1 , ^2 - 1 ] ■ The other parts of the Lemma follows from our construction. ■ 

Now, it is easy to prove Theorem|9]for r - \ using LemmafTOl Moreover, we can 
have a constructive algorithm to derive Pf from a given column p, which mainly has 
three steps: initialization, support reduction, and finalization. The step of initialization 
reduces the symmetry pattern of p to a standard form. The step of support reduc- 
tion is the main body of the algorithm, producing a sequence of elementary matrices 
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Ai , . . . , Ay that reduce the length of the coefficient support of p to 0. The step of final- 
ization generates the desired matrix P^, as in Theorem|9l More precisely, our algorithm 
written in the form of pseudo-code is as follows: 

Algorithm 1. Input p, which is a \ X s vector of Laurent polynomial with symmetry 
satisfying pp* — 1. 

1. Initialization: Lef q :- pU^p. Then Sq - [l.?,, -l.s,,z"'l.53, -z^'l^J, where all 
nonnegative integers ri, . . . , r4, si, . . . , S4 are uniquely determined by Sp. 

2. Support Reduction: Let Po := and J := 1. 

while (lcoeffsupp(q)l > 0) do 
if |coefFsupp(q)| > 1 then 

q := qBq, Ay := Bq, andPj := A}, 
else 

qis of the form: q = [0, 0, gi, -gajz"' + [£3, 0, gi, 92]; Lef 



By: = 











g;(i+z) 




gtd-r) 







2||gill 


2||gill 




g;(i-z) 
21lg,ll 





gju+z) 
2||gill 


G2 



where = UqiJ - 1,2. By is paraunitary, |coeffsupp(qBy)| = 0, and 

'>S(qBy) = [l,,,-l,„l,z-il,3_i,-l,-z-il,,_i]. 
Let E be a permutation matrix such that 

>S(qBy)£ = 
Ler Ay := EjE, Py := A}, andc\ = qAy. 
end if 
/ := / + 1 
end while 

3. Finalization: q - [£, 0, 0, 0] /or iome 1 x s'^ constant vectors £ in F. Let U := 
diag(f/f, Define Py := f/^. 

Output a desired matrix satisfying all the properties in Theorem\9\for r — I. 

4. Symmetric Complex Tight d-frame via Matrix Extension 

In this section, we shall discuss the application of our results on matrix extension 
with symmetry to d-band symmetric paraunitary filter banks in electronic engineering 
and to orthonormal multiwavelets with symmetry in wavelet analysis. In order to do 
so, let us introduce some definitions first. 

Let F be a subfield of C such that (I3.6l l holds. Let <p he a compactly supported 
d-refinable function in L2(K) associated with a low-pass filter aq : i-» F. Recall that 
its symbol is defined to be a()(z) := Y^kez «o(^)z*, which is a Laurent polynomials with 
coefficients in F. We denote its d-band subsymbols by ao;y(z) '■- Vd2/t£Z'3'o(T + dA:)z*^, 
y = 0, . . . , d - 1 . To construct a tight d-wavelet frame from (p, one has to design high- 
pass filters fl fl^ : Z t-^ F''^'^ such that jl-li holds. Let P{z) be the polyphase 
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matrix defined by 



ao;o(z) 
ai;o(z) 



ai;d-i(z) 



P{z) 

aL;o(z) ■■■ aL;d-i(z) 
Then it is easily seen that il.7\ is equivalent to 

nzYnz) = /d 



(4.1) 



(4.2) 



where each a,„;y is a subsymbol of a"' for m - 1, . . . , L; y = 0, . . . , d - 1, respectively. 

Symmetry of the filters is a very much desirable property in many appUcations. We 
say that the low-pass filter ao (or ao) has symmetry if 



a()(z) = z' 



(d-l)co 



ao(l/z) 



(4.3) 



for some cq e M such that (d - l)co e Z. If ao has symmetry as in ( I4.3l l and if 1 
is a simple eigenvalue of ao(l), then it is well known that the d-refinable function 
associated with the low-pass filter ao has symmetry; (p(co - ■) = <p. 
Under the symmetry condition in (14.3b . 



ao;y(z) = /''ao;e,(z-') r = 0, . . . ,d - 1, 
where y, 2y € F := {0, . . . , d - 1) and Ry, Qy are uniquely determined by 

(d - l)co -j^dRy + Qy with Ry eZ, QyE Y. 
Now, we can easily deduce symmetric Laurent polynomial via ( 14.4b . In fact, let 



(4.4) 



(4.5) 



bo;y(z) - 



ao;rk), 



y^Qy 



-i=(ao;y(z) + Z'^aQ;Q^{z)\ J<Qy 

^ (ao;y(z) - z*' ao;Q, (z)), y> Qy 



(4.6) 



where y - 0, . . . , d - 1 and each Ky is an integer so that the length of the coefficient 
support of a();y(z) + z'^''ao;() (z) is minimal. Then, one can show that each bo;y(z) is of a 
Laurent polynomial with symmetry. Note that the transform matrix U with respect to 
(14.61 ) is paraunitary. 

Next, we show that we can construct a vector of Laurent polynomial with symmetry 
from a low-pass filter a for the complex pseudo spline of type 1 with order (m, n) to 
which Algorithm[T]is applicable. 

Fixed m, « e N such that 1 ^ 2« - 1 ^ m. For simplicity, let oo - ci and ao 
be its symbol. From ( 14.61 ). we can construct a 1 x d vector of Laurent polynomial 
p(z) = [bo;o(z), . . . , bo;d-i(z)] from ao. However, pp* 1 when n < m. To apply our 
matrix extension algorithm, we need to append extra entries to p. It is easy to show that 



d-1 d-1 

1 - ^o(^ + 2jnl6f = 1 - ao;y(z=')a;^(z'), z = e 

j=Q y=0 
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where ao;y, 7 = 0, . . . , d - 1 are the subsymbols of oq. By Corollary |5] we have 



for some 27r-periodic trigonometric function bo with real coefficients. Hence, we can 
construct a Laurent polynomial ao;d(z) from bo such that ao;d(e"'^) = ^o(f)- Then, the 
vector of Laurent polynomials q(z) - [ao;o(z), ■ • ■ , ao;d-i (z)> stoAz)] satisfies qq* - 1 . 

For m = 2n, by Corollary |5l |^o(^)P = C2„,2n-i[sin^(^/2)/d^]^"-\ W ao;d(z) can be 
constructed explicitly as follows: 



Then ao;d(z) is symmetric, let bo;ci(z) = ao;d(z)- Then p := [bo;o, ■ ■ ■ , bo;d] is a 1 x (d + 1) 
vector of Laurent polynomials with symmetry satisfying pp* = 1 . 

Ifmi^ 2n, ao;diz) is a not necessary symmetric, we can further let bo;d(z) - (a();d(-z)+ 
ao;d(l/2))/2andb();ci+i(z) = (a();d(z)-ao;d(l/z))/2. In this way, p := [bo;o, ■ ■ ■ , bo;d, bo;d+i] 
is a 1 X (d + 2) vector of Laurent polynomials with symmetry satisfying pp* = 1 . 

Consequently, we can summerize the above result as follows: 

Theorem 11. Let m,n e N be such that I ^ n ^ m. Let ao be the low-pass filter 
for the complex pseudo spline of order (m, n) defined in (11.121 1 and ao(z) be its symbol. 
Then, one can derive Laurent polynomials ao;cj(z), . . . , ao;L(z), d ^ L ^ d + 1 such 
that pao := [ao;o, ■ ■ ■ , ao;d-i, ■ ■ ■ , ao-iiz)] satisfies PaoPao = ^' "^here ao-o, ao-a-i are 
subsymbols of ao- Moreover, one can construct a{L + l)x(L+ 1) paraunitary matrix 
U such that PaoU is a vector of Laurent polynomial with symmetry. In particular, for 

m = 2n, a();d(z) = ^/C2n,2n-\ { ^ ij^^ ) ^ and L - d. 

Now, applying Theorem|9] we have the following algorithm to construct high-pass 
filters ai, . . . , aL from a low -pass filter ao for a complex pseudo spline of type I with 
order (m, «). 

Algorithm 2. Input low-pass filter ao for a complex pseudo spline of type I with order 
(m, n), I ^ 2n — I ^ m. ao satisfies ( 14.31 ). 

(1) Construct Pao(z) and U as in Theorem 1771 such that p := PaoU is a I x (L + 1) 
row vector of Laurent polynomials with symmetry (L — d when m — 2n while 
L = d -H 1 when m + 2n. 

(2) Derive Pj. with all the properties as in Theorem\9[for the case r — I from p by 
Algorithm\l\ 

(3) Letf :— PcU* —: (am;y)osjm^L,o^y^d-i OS in ( 14. 11 1. Define high-pass filters 



d-1 



l-2fe(^ + 2>/d)|' = |feo(da'. 





d-1 



m - \, . . . ,L. 



(4.7) 



Note that we only need the first d columns o/PgU*. 
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Output a symmetric filter bank {ao, ai , . . . , 3ll} with the perfect reconstruction property, 
i.e. VizyPiz) — Id. All filters a™, ot = 1, . . . , d - 1, have symmetry: 

a,„(z) = fi,„z^^"'-^°a„,(l/z), (4.8) 

where c'" :— km + co e M and all s„, e {-1, 1), k,„ e Z/or m — I, . . . ,L are determined 
by the symmetry pattern of Pg as follows: 

[l,siz'",...,SLz'''VSp:^SPe. (4.9) 

Proof. Rewrite P^, = (b,„;y)o^m,y^L- Since has compatible symmetry as in (|4.9l l. we 
have Sbm-y - e„,z* >Sb();y for m = 1, . . . , d - 1. By ( I4.6l i, we have Sbo-y = sgn(gy - 
^-j^Ry+ky^ y = 0, . . . ,d - 1, where sgn(x) - 1 for x ^ and sgn(x) = -1 for x < 0. 
Then, we have 

Sb„r,y = sgniQy - r)e„,z'''^^'^*"', (4. 10) 

By (14.101 1 and the transformation matrix U* with respect to (I4.6l l. we deduce that 



(4.11) 



This implies that >Sa,„ - e,„z'^**"'^^"^ which is equivalent to ( I4.8l l with c,„ :- k„, + cq 
for m - 1 , . . . , L. 

Since the high-pass filters a\,. . .,3.1 satisfy (14. 8t . it is easy to verify that , . . . 
defined in ( |1.4t also has the following symmetry: 

iA'(ci--) = eiiA', iA^(c2 - ■) = £2'A^ ip^{cL--)^SL>p\. 

In fact, by (14.8b . (11.4b and the symmetry of (p, we have 

r(-cl^) = a„,(e'0?(-^) = fi™e-'<^^'"-''°'a„,(e-'f)e-''^»f^(^) 



(4.12) 



„-'Xdc„,)fE 



a;„(e '^)0(^) = e„,e 



-'■(df„,)f 



iA'«(d^). 



This impUes - s„,e^"'"'^iff'"{^), which is equivalent to ( 14. 12b . 

In the following, let us present some examples to demonstrate our results and illus- 
trate our algorithms. 

Example 1. Consider dilation factor d = 2. Let m = 4, « = 2. Then ^43(3') = 1 -1- 4^ -H 
4y^. The low-pass filter a^'^ with its symbol ao for the complex pseudo spline of order 
(4, 2) is given by 

\4 r/ 



ao(z) 



1 +z 



V6 



- +12 + 



V6 



1 



V6 



1 



Note that coeffsupp(a()) = [-3,3] and a(z) = a(z '). In this case, m = 2n, by Theo- 
remim we can obtain Pa„ - [ao;o(z), a{);i(z), stoaiz)] as follows: 

nr. 

ao;o(z) = - ^ (3z + 8 V6/ - 6 + 3z"'); 
ao;i(z) = + 8 V6/ - 26 + 5z-')(l + z"'); 



a0;2(z) 



160 

^(z-2 + z-')(l-z-'). 
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Pao is already a 1 x 3 vector of Laurent polynomials with symmetry, i.e., Sp^o — 
[l,z"', -z"']. Applying Algorithm|2] we can obtain its 3 x 3 extension matrix P* - 
[p*^,p;j,p;j with Pa, := [ai;(),ai;i,ai;2] and pa^ := [a2;o, a2;i, a2;2] as follows: 

a.o(z)^ -^^^-^^'^ 3.-8V6,.42.3z-); 

au(z) = -^l^f "'^ (5z - 8 V6/ - 74 . 5z-') (1 . z^); 

ai;2(z) = "^\f^;^'^ (5z - 16 V6/ - 58 + Sz-^ d " ^-'); 

JO 

a,,(z) = (5z - 4 V6/ - 22 + Sz"') (1 - z^); 

a2;2(z) = (5z - 12 V6/ - 26 + Sz"') (1 + z"'). 

Note that the symmetry of P^. satisfies SPe - [1,1, -l]''^>Spa„ and the coefficient sup- 
port of Pg satisfies coeffsupp([Pj,]:,y) c coeffsupp([pa„]j) for j = 1,2,3. The high-pass 
filters ai, a2 constructed from Pa, and Pa, via jA.lj are then given by: 

a,(z) . - ^^^^^^;^'\ 5z . 26 . 2 V6/ . 5z-')(z - 2 . z^)^; 

a2(z) = ^(^ + ^'^ 5^ + 16 -H 2 V6/ -H 5z"')(z - 2 -h z"')(z - z"'). 
560 

We have ai(z) = ai(z"') and a2(z) - -a2(z"'). Let be the 2-refinable function 
associated with the low-pass filter ao. Let i/'',!^^ be the wavelet functions associated 
with the high-pass filters ai, a2 respectively by ( 11.41 ). Then 4>(--) = (p, (/''(-■) - 4'^ and 
= -i/'^. The graphs of (p, i/'^ are as follows: 



Figure 4.1: The graphs of (p,i// ,if/ (left to right). Real pai1:solid Hne. Imaginary partidashed line. 

Example 2. Consider dilation factor d = 3. Let m = 4, n = 2. Then ^4,3(3') =1-1- 
+ 64^^). The low-pass filter a^'^ with its symbol ao for the complex pseudo spline 
of order (4, 2) is given by 



/ 1 



ao(z) 



-H 1 -HZ 



( 4 2V5 



1 (n 4V5. 



+ 

3 3 



( 4 2V5 
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Note that coeffsupp(a()) = [-5,5] and a(z) = a(z ^). In this case, m - 2n, by Theo- 
remfm we can obtain Pa„ - [ao;o(z), ao;i(z), ao;2(z), ao;3(z)] as follows: 

a,);o(z) = (lOz + 27 V5/ - 20 + lOz"'); 

ao;i(z) = ^(-(4 + 2 VsOz"^ + 30z-' + 60 + 6 Vs/ - (5 + 4 VsOz); 

Pi 

ao;2(z) = ^(-(5 + 4 V5/)z-2 + (60 + 6 VsOz"' + 30 - (4 + 2 V5/)z) 

2 VTo , 

ao;3(z) = — ^(z-2 + z-')(l-z). 

Note that ao;i(z) - z"'ao;2(z"')- Let p - Pa„U with U being the paraunitary matrix 
given by 

10 





U : = 



^ ^ 
J- -J- 











Then p is a 1 x 4 vector of Laurent polynomials with symmetry pattern satisfying Sp 



[hz- 



']. Applying Algorithm|2] we can obtain a 4x4 extension matrix P* 



[Pao'Pa.'Paj'Pas] with pa, := [ai;o, ai;i, ai;2, ai;3], Pa2 [a2;o,a2;i,a2;2,a2;3], and 
Pa3 ■- [a3;o, a3;i, a3;2, a3;3]. The coefficient support of Pf satisfies coeffsupp([Pf]:j) c 
coefFsupp([pa„] /) for j = 1,2,3,4. The high-pass filters ai , a2, a3 constructed from paj 
, Pa,, and Pa, via ( I4.7l i are then given by: 



V19178 , 
^'^^^ = 4660254^^'^^^^^'^^^^^ 
V2 18094 , 
i /66j6i4 



a,(z)-^^^{h(z)-z'b,{z-')). 



where 



biiz) = (-172 - 86 / V5)z5 + (-215 - 172 / V5)z'* - 258 / VSz^ 

+ (1470 + 1224 / Vs) + (i860 + 2328 / Vs) z - 3036 / Vs - 2943 

b2(z) = (-652 - 326 / Vs) z^ + (-8 15 - 652 / Vs) z^ - 978 / VSz^ 
+ (l832 ;• V5 + 1750)z^ + (3508 / Vs + 3020)z 

bi(z) = (4 V5 + 10/)z^ + (5 V5+20/)z'* + 30iz^ +(-53 V5-260i)z^ 

We have ai(z) = ai(z"'), a2(z) = -a2(z"'), and a3(z) - -z-'a3(z"'). Let (p be the 3- 
refinable function associated with the low-pass filter ao- Let iff\i/j', if/^ be the wavelet 
functions associated with the high-pass filters ai,a2,a3 respectively by ( 11.41 ). Then 
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<^(-) = 0, iA'(-) = >P\ >P\-) 

graphs of (p, i/'', (A^, and . 



-ijp', and - ■) - -ijr' . See Figure 2 for the 



Example 3. Consider dilation factor d - 3. Letm -5,n-2. Thenf ^^Cy) = 1 + yy + 
^y^). The low-pass filter with its symbol a^' ^ for the complex pseudo spline of 
order (5, 2) is given by 



ao(z) = 



+ 1 +z\ 



In this case, m = 2« - 1, we have 



2 



V30 ^ 
1 



1 



1 

- + 

z 



(13 2V30, 



-i\z 



1 - 2^ ao;,(z)*ao;^(z) = ——(1 - z)(l - z-')(4z^ + 106z^ - 54U 



y=0 



19683 

+ 2320 - 541z-' + 106z"2 + Az"^), 



By Riesz lemma, we can factorize the above positive trigonometric function and obtain 
ao;d(2) = co(z-ri)(z2-r2Z+r3)(l-z) withri = ci+ - 1, = C2 + q andrs = |c2p, 
where ci e M, C2, C2 e C are the roots for the polynomial + 2123?^ - 553>y + 1054 
and CO = .,Jf e R. 

243 ^/-rirj. 

Letpa„ = [ao;o(z),ao;i(z),ao;2(z),z"^-^ao;3(z),z"^^a();3(z"')]. Then paoPa„ = 1. 

Note that a();i(z) - z"^ao;2(z"')- Let p - PaoU with U being the paraunitary matrix 
given by 



1 1 

V2 V2 
1 l_ 

V2 V2 





u 















1 

f 








1 



V2 V2 

Then p is a 1 x 5 vector of Laurent polynomials with symmetry pattern satisfying 
Sp - [l,z"',-z"', 1,-1]. Applying Algorithm|2] we can obtain a 5 x 5 extension 
matrix P^. The coefficient support of satisfies coeffsupp([Pe]:j) c coeffsupp([paJ;) 
for j = 1 , . . . , 5. From P^., we can derive high-pass filters ai , . . . , a4 via (I4.7l i as follows: 

ai(z) = bi(z) + bi(z-'y, a2(z) = biiz) - biiz-'); 
aiiz) = bi(z) + z'biiz-'y, a4(z) = b^iz) - z^b^iz-'); 



21 



where 

biiz) * -0.38736 + 0.26298/ - (0.00307 + 0.002040z^ - (0.00744 + 0.00660/) 

- (0.00951 + 0.01459/)z'* + (0.01372 - 0.02716/) 

+ (0.06342 - 0.03967/) + (0.13657 - 0.04144/) z; 
b2(z) ~ (0.00601 + 0.00672/) + (0.01353 + 0.02017/) + (0.01353 + 0.04034/)z'* 
-(0.01370- 0.03024/) z"* 

- (0.07035 + 0.0001 l/)z^ - (0.14727 + 0.03064/) z; 

bi(z) ~ - (0.01063 + 0.01 189/) z^ - (0.02392 + 0.03567/) z^ - (0.02392 + 0.07133/)z'* 

+ (0.05891 - 0.08433/) z^ - (0.00043 - 0.20320/) z^; 
b4(z) ~ (0.00417 + 0.00467/) z'' + (0.00939 + 0.01400/) z^ + (0.00939 + 0.02800/)z'* 

- (0.02126- 0.02655/) z^ - (0.23403 + 0.31 140/) z^ 

The coefficients of bi{z), . . . , b4{z) are rounded from exact explicit solutions (too long 
to be presented here). We have ai(z) - ai(z"'), a2(z) = -a2(l/z), a3(z) = a3(l/z), 
and a4(z) = -z^a4(z"'). Let <p be the 3-refinable function associated with the low- 
pass filter ao. Let ij/^ , . . . , iff* be the wavelet functions associated with the high-pass 
filters ai, . . . , a4 respectively by ( fL4] i. Then (f>{—) - (p, - i^', iA^(— ) = -ifr^, 

- ■) = ip^^ and - ■) = -lA^ • See Figure 3 for the graphs of 0, if>^,i//^, ifr" and 






Figure 4.3: The graphs of (fi,tff' (left to right, top to bottom). Real part:solid line. Imaginary 

partidashed line. 



Example 4. Consider dilation factor d = 3. Let m = 4, « = 2, and the low -pass filter 
as in Example 2 with its symbol ao for the complex pseudo spline of order (4, 2) 
is given by 



/ 1 



ao(z) 



-H 1 -HZ 



( 4 2V5 
3 3 'j 



1 (n 4V5. 



H 1 

z \3 3 



( 4 2V5 
3 3 'j 
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By Theorem [8] the shifts of 4> are linearly independent. Hence, there exist compactly 

supported dual refinable functions in e ^2(1^) for 0, i.e., 0(d^) = flo(^)0(^) for some 
low -pass filter flo and {<p, 0(- - k)) = 6k,ke 7L. Here, we provide a low-pass filter aofe) 
for 4> as follows: 



where 



ao(z) = 



329387 209689 / 102661 5464379 
b(z) = + iV5-| -I- W5|z 



2754 1377 
/ 177727 551620 



\ 2754 4131 



\ 816 22032 

/2967467 1034833 



+ - 



375253 158555 



4131 
24103 



15147 



V5 



z + 



\ 22032 22032 
/ 24620753 29059 
\ 727056 " 22032 



V5 



3366 1 16 727056 



See Figure 4 for the graph of the 3-refinable function (p associated with the low-pass 
filter aofe)- Note that = 0(— )■ 





Figure 4.4: The graph of ^. Real part:left. Imaginary part:right. 
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